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Abstract 

We compute numerically the yrast line for harmonically trapped boson sys- 
tems with a weak repulsive contact interaction, studying the transition to a 
vortex state as the angular momentum L increases and approaches A^, the 
number of bosons. The L = N eigenstate is indeed dominated by particles 
with unit angular momentum, but the state has other significant components 
beyond the pure vortex configuration. There is a smooth crossover between 
low and high L with no indication of a quantum phase transition. Most strik- 
ingly, the energy and wave function appear to be analytical functions of L 
over the entire range 2 < L < N. We confirm the structure of low-L states 
proposed by Mottelson, as mainly single-particle excitations with two or three 

units of angular momentum. 

PACS numbers: 03.75.Fi, 05.30.Jp, 67.40.Db 
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The low-lying excitations of atomic Bose Einstein condensates in harmonic traps 
are of considerable experimental and theoretical interest Recently, Mottelson proposed 
a theory for the yrast line of weakly interacting A^-boson systems 0, i.e. the ground states 
at nonvanishing angular momentum L. Physical arguments led him to assume that the 
yrast states are excited upon acting on the ground state |0) of vanishing angular momentum 
with a collective operator Q\ = Y.p=i that is a sum of single-particle operators acting 
on the coordinates Zp = Xp + iyp of the p^^ particle. For angular momenta L <^ N the 
yrast states are found to be dominated by quadrupole (A = 2) and octupole (A = 3) modes. 
Assuming a vortex structure of the yrast states with L ^ N then led to the prediction of 
a quantum phase transition in Fock space when passing from the low angular momentum 
regime L -C to the regime of high angular momenta L ^ N. The reason for this is 
behavior is the approximate orthogonality of the collective states Q\\0) and the single- 
particle oscillator states of the vortex line in the regime N^^'^ <^ L. These results have been 
obtained for harmonically trapped bosons with a weak repulsive contact interaction. The 
case of an attractive interaction has been studied by Wilkin et al. [Q. In this case the total 
angular momentum is carried by the center of mass motion, and there are no excitations 
corresponding to relative motion. This is not unexpected since internal excitations would 
increase the energy of the yrast state. 

It is the purpose of this letter to present an independent numerical computation of the 
yrast line and to compare with Mottelson's results ||^. In particular we want to focus on 
the transition from low to high angular momentum yrast states. The investigation of this 
transition is of interest not only for the physics of Bose-Einstein condensates. Localization 
in Fock space is under investigation also in molecular and condensed matter physics . 
The numerical computation has the advantage that it does not rely on the assumptions 
made in the analytical calculation. However, with our numerical methods it is limited to 
angular momenta below about L ^ 50. Most interestingly, our numerical results suggest 
that the yrast line and the corresponding wave functions can be presented by rather simple 
analytical expressions. 

2 



Let us consider bosons in a two-dimensional harmonic trap interacting via a contact 
interaction^. We are interested in the yrast hne in the perturbative regime of weak interac- 
tions. Note however that experimental studies of trapped condensates are often in a regime 
where the interaction energy is comparable to the trapping potential, and this may introduce 
qualitatively different physics. We write the Hamiltonian as 

H = Ho + V. (1) 

Here 

Ho = huj^ia]aj (2) 

3 

is the one-body oscillator Hamiltonian and 

^ = 9 Vijkia\a]akai. (3) 

the two-body interaction. The operators dm and annihilate and create one boson in the 
single-particle oscillator state \m) with energy mhu and angular momentum mh, respectively 
and fulfill bosonic commutation rules. The ground state energy is set to zero. Up to some 
irrelevant overall constant the matrix elements are given by Vijki = 2^''^''{k + {iljlkW.y^'^ 
and vanish for i + j ^ k + I. For total angular momentum L the Fock space is spanned 
by states \a) = \no,ni, . . . ,nk) with Ei=o,fe'^i = a]aj\no,ni, . . . ,nk) = nj\no,ni, . . . ,nk) 
and J2j=o,kj^j = Here rij denotes the occupation of the j^^ single particle state 
For vanishing coupling g the basis states are degenerate in energy, and the problem thus 
consists in diagonalizing the two-body interaction V inside the Fock space basis. To set up 
the matrix we act with the operator (^) onto one initial basis state with angular momentum 



L and onto all states created by this procedure until the the Fock space is exhausted [10 



The resulting matrix is sparse, and the yrast state is computed using a Lanczos algorithm 



Tl| . We restrict ourselves to L < 50 corresponding to a maximal Fock space dimension of 
about di ^ 2 ■ 10^ 



^The results obtained below extend to the three dimensional problem for L = 
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The yrast line, i.e. the ground state energies as a function of the angular momentum 
may be written as E{L) = Lhuj + gei. Fig. |1] shows the L-dependence of the energies for 
systems of = 25 and 50 bosons. The energies simply decrease linearly with increasing 
angular momentum for L < N . In fact to machine precision, the energy function is found 
to described by an algebraic expression. 

At fixed angular momentum L and ioi L <^ N the energies gei increase as expected with the 
square of the number of bosons N . Notice in the figure that there is a kink in the slope at 
N = L. This is a hint of condensation into a vortex state: in macroscopic superfiuids, the 
state for L = N would have a condensate of unit angular momentum and would be lower in 
energy than neighboring yrast states. 

We next investigate the structure of the wave functions of yrast states. We would like 
to know how complex the states are and how well they can be described by single-particle 
operators acting on simple states. To address the question of the complexity of the states 
in the Fock basis, we take the wave function amplitudes c^^^ in the Fock representation of 
the state 

and compute the inverse participation ratio ||12| 



dL 
a=l 

The II is the first nontrivial moment of the distribution of wave function intensities jc^^-'p. 
Its inverse measures the number of basis states \a) that have significant overlap with the 
yrast state \L). Fig ^ shows a plot of 1//^ and the Fock space dimension as a function of 
angular momentum L for a system of = 50 bosons. The is seen to be much smaller 
than the dimensionality of the Fock space. Even where the participating is greatest, at 
midvalues of L, only about 30 states are active participants. A similar behavior of quantum 



non-ergodicity has been found previously in numerical studies [10|. Notice that the inverse 



participation ratio decreases strongly as = L is approached. This shows that the yrast 
state becomes simpler, again hinting at the formation of a vortex condensate. Examining the 
coefficients for = 25 in detail, the largest amplitude at L = is in fact the vortex state, 
la) = |OA^O...O >, but it has less than half the probability of the complete wave function. 
Interestingly, our numerically obtained yrast state \L = N = 25) agrees with the conjecture 
given by Wilkin et al. i.e. |L = A^) = nJLi(^p - Zc) |0) with = N'^ J2p=i Zp being the 
center of mass. Based on our numerical wave functions, we can generalize this conjecture. 
We believe that all the yrast states for 2 < L < A^ are given by the formula 

1^) = E izp,-Zc){zp^-Zc)...izp^-Zc)\0> (5) 

Pl<P2<---<PL 

We have verified that this formula is correct (up to machine precision) by comparison with 
the numerically obtained yrast states for A^ = 25. Since the operator acting on the ground 
state is translationally invariant no quanta of the center of mass motion are excited. Notice 
that there is a natural termination of the construction at L = A^. 

To further examine the structure of the yrast states \L) we show a plot of the occupation 
numbers n^^^ = (L|djaj|L) for j = 0, 1, 2, 3 in Fig. ^for a system of A^ = 50 bosons. At very 
low angular momenta the yrast states are dominated by single particle oscillator states with 
two or three units of angular momentum. This is in agreement with Mottelson's results P|. 
However, at larger angular momentum L the dominant fraction is carried by single particle 
states with one unit of angular momentum. Note that the occupation numbers n^^^ are 
very small for j > 3. This analysis confirms the results found for the inverse participation 
ratio. Note also that the observables n^-^^ are very smooth functions of L. If there were a 
quantum phase transition at large L ^ A^/2, we would expect to see some precursor in these 
observables. 

In conclusion, our numerical study strongly indicates that there is no quantum phase 
transition to a vortex state for trapped condensates in the limit that the interaction poten- 
tial is small compared to the oscillator frequency. The strongest evidence is the apparent 
existence of analytic expressions for the energies and the wave functions on the yrast line for 
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2 < L < N One might speculate that these states are contained in a dynamical symmetry 
group, but we have no idea how this might come aboutQ. We have also examined the struc- 
ture of the yrast states and the matrix elements between them, finding that the observables 
vary smoothly with L, for L not too small. 

We acknowledge conversations with B. Mottelson. This work was supported by the Dept. 
of Energy under Grant DE-FG-06-90ER40561. 



^We note also that there is another symmetry group |T^, SO{2,l) , that produces relationships 
between energies of different states within a single L subspace. 
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FIGURES 




FIG. 2. Inverse participation ratio (circles) and Fock space dimension (11 (squares) as a 
function of angular momentum L for a system of = 50 bosons. 
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FIG. 3. Occupation numbers of the lowest single particle oscillator states as a function of 
angular momentum L for a system of = 50 bosons. (j = 0: diamonds, j = 1: squares, j = 2 
circles, j = 3 triangles) 
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